Fluctuation-induced forces between inclusions in a fluid membrane under tension 
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We develop an exact method to calculate thermal Casimir forces between inclusions of arbitrary 
shapes and separation, embedded in a fluid membrane whose fluctuations are governed by the 
combined action of surface tension, bending modulus, and Gaussian rigidity. Each object's shape 
and mechanical properties enter only through a characteristic matrix, a static analog of the scattering 
matrix. We calculate the Casimir interaction between two elastic disks embedded in a membrane. 
In particular, we find that at short separations the interaction is strong and independent of surface 
tension. 
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Biological membranes are formed by bilayers of am- 
phiphilic lipid molecules with embedded macromolecules 
such as proteins jlj]. While inclusion can dramatically 
modify the flexibility, surface tension, and spontaneous 
curvature of biomembranes [^) , the interactions between 
embedded proteins mediated by the membrane, in turn, 
depend on its mechanical properties J^-Q]. Generally, 
these interactions depend on the separation between 
the inclusions and can be categorized as elastic and 
fluctuation- induced Q . The former is due to static elastic 
deformations induced in the membrane, while the latter, 
the Casimir-like force, is due to the modified entropy of 
the membrane fluctuations. 

Despite being the subject of intense research for al- 
most two decades^, ^, ||, fluctuation- induced forces 
are just beginning to be understood. Previous theoreti- 
cal research only considered fluctuations in the separate 
presence of either surface tension or the bending rigidi- 
ties, but not the two combined |To| , as occurs in real- 
ity. While surface tension is negligible for free-floating 
membranes it is finite for a membrane enclosing a cell 
with excess osmotic pressure, or for membranes under 
tension in general|lT|. In such cases, the surface tension 
cr , as well as the bending (k ) and Gaussian (~k = hqKq) 
rigidities of the membrane are all non-zero. For such 
membranes, the static elastic interaction between inclu- 
sions is exponentially cut-off at distances larger than the 
characteristic length ||, [l^] 

£ =a^ = (Wo) 1/2 - (1) 

In this letter, we present an exact method that allows 
calculations of membrane-mediated fluctuation-induced 
forces between any number of arbitrarily placed elastic 
inclusions of any shape for finite £q. It is based on the 
technique previously employed to calculate the Casimir 
forces in antiferromagnets fl3f . It can also be regarded as 
a version of the scattering-matrix approach j^j, general- 
ized to entropic forces for Hamiltonians with quartic in 
derivatives terms representing bending energy. We find 
that the Casimir energy between objects embedded in 
membrane could be expressed in terms of the response 
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FIG. 1. (color online) The Casimir free-energy (|l2j) scaled 
with fourth power of the distance as a function of R/a for 
k/ko = 10 4 (symbols) and k/ko — 10 — 1 (solid lines), with 
aoa = a/lo as indicated. We set «o = — kq, k = — K. Dashed 
black lines indicate asymptotic large- R dependences evaluated 
with Eq. (|l6|). Plots are horizontal (AT <x R~ A ) only for 
R < £q = ctg 1 ; for larger R/lo the Casimir energy decays 
as AT oc R~ 8 , see Eq. (^). Inset: The Casimir energy of 
two disks in the regime dominated by the bending energy 
(to/a = 10 4 , k/k = 10 4 ) divided by the PFA estimate @. 
The quantity m indicates the multipole order of truncation. 

of individual objects to the fluctuating field. An object's 
shape and material properties enter only through the co- 
efficient matrix representing the object's response, which 
is a static analog of the scattering matrix. 

In particular, we compute fluctuation-induced interac- 
tions between two elastic disks. We specifically exam- 
ine the unexplored parameter range where £q is between 
the disk radius a and the center-to-center separation dis- 
tance R, a < £o < R, so that both the surface tension and 
bending energies are relevant [Fig. . In the limit of very 
short separations, we find a very good match between our 
results and those obtained through Derjaguin|l^] or prox- 
imity force approximation (PFA) [sec the inset in Fig. Q . 
In the large separation limit, R 3> a, we derive explicit 
asymptotics for the Casimir interaction, and check them 
against the known results ||, ||, |[ ^| for cases dominated 
by either surface tension or bending rigidities. 

For a non-zero surface tension, we find that the Casimir 
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energy becomes strongly suppressed at distances R larger 
than the characteristic length £g, Eq. ([j]). Note that the 
Casimir energy retains the power-law asymptotic form 

/3F C = -Aa n /R n , p = l/k B T, R > a, (2) 

with the exponent increasing from n — 4 at R < £q 
(bending-energy-dominated regime) to n — 8 at R > £q 
(tension-dominated regime). This is illustrated in Fig. [j] 
for two inclusion stiffnesses, with a number of differ- 
ent surface tensions. Depending on the parameters, the 
Casimir free energy scaled with the fourth power of the 
distance has either constant or oc 1/i? 4 asymptotics. 

As seen in the figure, all curves with the same ratio of 
k/kq merge at short separations, while at larger separa- 

— 1/2 

tions the value of lo /a oc tr plays an important role. 
We also find a surprising effect corresponding to the dif- 
ference between the Gaussian rigidities of the membrane 
and the inclusions. In the absence of line tension 
for the inclusion boundaries [a — a a in Eq. (^)], the 
Casimir force at large distances becomes exponentially 
small when the Gaussian rigidities of the inclusions and 
the membrane coincide, re = reo, regardless of re. 

We write the energy of a given configuration as a com- 
bination of that of the inclusions (thin isotropic elastic 
solids characterized by in-plane Lame coefficients and 
bending rigidities) and the fluid membrane outside the 
inclusions (defined by the sum of the surface tension and 
bending energies). The energy is expanded to quadratic 
order in the displacements with respect to the equilib- 
rium membrane configuration assumed to be planar in 
the z — plane. The result for the membrane with in- 
clusions has the standard Helfrich form |l7]] , 



U= j d 2 r|(V U ) 2 + f(V 2 U ) 2 + «[< K < a -« y ) 2 ], (3) 



where primes (') denote the partial derivatives with re- 
spect to x or y as indicated, and the integration is done 
over the total projected area A. The coefficients k and 
k in Eq. (|3|) are position-dependent, e.g., k = «(r), 
r = (x,y). Thermodynamical stability requires that 
a, re > and -2re < re < in Eq. (§). 

The first term in Eq. (||) has the standard form of a 
surface-tension contribution. However, inside inclusions, 
it represents the elastic energy associated with the in- 
plane stress induced by the membrane surface tension. 
In the absence of line tension at the inclusion boundaries, 
the diagonal components of the equilibrium stress tensor 
in the inclusions coincide with the surface tension of the 
membrane, which gives a = <tq in Eq. (3). 

The terms with re and re in Eq. (3) represent the 
bending energy contributions associated with the mean 
(re = kq + Arei(r)) and Gaussian (re = reo + Arei(r)) cur- 
vatures respectively. While reo and reo correspond to an 
unperturbed membrane, rei(r) and rei(r) are position- 
independent inside, and vanish outside of the inclusions. 



The partition function Z of an inhomogeneous mem- 
brane can be found as a Boltzmann sum over all mem- 
brane configurations u(r). With quadratic Eq. (||), the 
free energy, T = -^InZ = (2/3)" 1 £„ ln(/3£ n ) + 
const, is expressed in terms of the eigenvalues E n of 
the hermitian "Hamiltonian" operator H = H T ob- 
tained as the second functional derivative of Eq. (^) over 
w(r). The corresponding set of orthonormal eigenfunc- 
tions it„ = it n (r), Hu n — E„u n , is complete in the space 
of the functions with continuous second derivatives on A. 
We define the Green's function (GF) of the operator T-L, 



G^G(r,r') = ^ 



u n (v)u n (r') 



(4) 



which obeys the usual equation H r G(r,r') = S(r — r'). 
For a uniform membrane, the term with re in the inte- 
grand of Eq. (|) is a total derivative; the correspond- 
ing energy operator is Ho = reo(V 4 — a 2 V 2 ), with a 2 , = 
co/reo The uniform- membrane GF is a combination of 
GFs for the Laplace and the Helmholtz equations |13] , 



Go(r) 



[K (a r) +ln(a r)], r 



(5) 



Note that the operator H = 7io + XV, the eigenfunc- 
tions u n = u n , the energies E n = E^ , and the free energy 
T = T\ depend on the coupling parameter A, with A = 0, 
corresponding to a non-perturbed (uniform) membrane. 
We use the Hellmann-Feynman theorem ]l8| to write the 
part of the free energy associated with the inclusions, 
AJ-\ = J-\ — J-q, as an integral over the parameter A, 
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dX' Tr(VG A <)- 



(6) 



For k > 1 membrane inclusions, we further decompose 
the free energy into the part associated with individual 
inclusions, X)f=i ^*^a > an( ^ ^ ne Casimir energy proper, 



Jb ee A^ A - V AT 7 - 



(0 



(7) 



To this end we write V — J2i=i Vl> where the operator V;, 
/ = 1, 2, . . . , fc, is only non-zero inside the corresponding 
inclusion. Then, if we expand the GF in the integrand of 
Eq. (^) in the perturbation series over A', the n-th term 
is a loop with the unperturbed Go connecting exactly n 
operators V/ . We group together subsequent terms corre- 
sponding to the same inclusion I by introducing the exact 
GF for a single inclusion, 



G 



(0 



Go — AGqV/Go + A G, : V;G,,V;G|, 



(8) 



Respectively, the Casimir free energy (Q) comprises the 
terms with at least two unequal inclusion indices p9|, 
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n>l 



(-AT 
2n 



{h} 



•V;„G 



(In) 



(9) 
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where the n-th term involves the summation over n in- 
clusion indices 1 < /, < k, with the neighboring indices 
different, li, l n l\. Physically, this can be in- 

terpreted as fluctuations' back-and-forth "hopping" be- 
tween the inclusions. To evaluate the resulting series, we 
introduce the k x k operator-valued matrix E with ele- 
ments S;j/ = (1 — 5/;/)V/G^' (no summation), and write 
the inner sum in the n th term of Eq. (^) as 



TrE™ 



y^TrE; li2 E; 2 ; 3 

{li} 



Z ln _ lln Vi nh ,n>l, (10) 



where neighboring indices are automatically different. 
Now, the Casimir free energy becomes 



1 



#Fo = -Trlog(l + AE) 



(11) 



This equation is our main general result: it is exact, re- 
mains finite in the limit A — > oc, and can be applied to 
calculate the Casimir forces between a finite number of 
compact objects of arbitrary shape and separation. The 
result can also be recast |ll| in a form similar to that of 
the scattering matrix approach employed to calculate the 
electromagnetic (EM) Casimir interaction jj3. 

We now focus on the case of two uniform circu- 
lar disks embedded in a membrane. The correspond- 
ing matrix E is 2x2, with no diagonal elements, 
and, therefore, the terms with odd powers of A in the 
expansion of Eq. (|ll| ) disappear. The matrices E", 
with even powers n — 2s, are diagonal, since E 2 — 

diag(ViG^ 1) V 2 Gi 2) ,V2Gi 2) ViGi 1) )- The two matrix el- 
ements give equal contributions to the trace, 



0?c 



Tr[A 2 ViG^V 2 G\ 2) ] s 



| Tr log(l - X A V l G^ ) \>2G^ > ). 



:( 2 h 



(12) 



For actual calculations, we construct the exact single- 
disk GFs Gf, 1 = 1,2, as a series in polar coordinates, 
including the terms corresponding to both the Laplace 
(oc r ±m ) and the Hclmholtz [oc K m (ar), I m (ar)\ equa- 
tions [cf. Eq. (^)], with the asimuthal quantum num- 
ber to < TO max . This requires four boundary conditions 
on the circumference: continuity of the function, normal 
derivative, as well as of the following two quantities, 



Q3 = ad r u — nd r (V 2 



-d r [-u 9t 



Q 4 = kV 2 u + - ( -u'g, 
r \r 



d r u 



(13) 
(14) 



Here u'g is the second derivative over the polar angle 
with respect to the center of the disk. The same bound- 
ary conditions are used to evaluate the matrix elements 
of the operators V; in Eq. (fl2| ). This gives the argument 
of the logarithm in Eq. (\12\) as a 2m nl&x x 2m max ma- 
trix (mode doubling corresponds to Helmholtz/Laplace 



components), which is non-diagonal since the GFs are 
expanded with respect to two different centers. 

At large separations, it is sufficient to keep the terms 
up to quadrupole (m max = 2). As the distance R between 
the disks decreases, higher order multipoles become rel- 
evant. Generally, when the distance between the edges 
of the disks, H = R — 2a, is small, a large number of 
multipoles, m max > a/H, is required for convergence. In 
this regime the Casimir energy becomes large, and inde- 
pendent of the surface tension (see Fig. |l|) . 

The short separation asymptotic form of the Casimir 
energy can also be evaluated within PFA|D|. Here, 
the interaction between curved edges is expressed as a 
sum over infinitesimal straight line segments approxi- 
mated as parallel. We found that the Casimir energy 
per unit length for two half-planes is (3Tc/L = f/H, 
with / = 7r/24 in the limit dominated by the surface ten- 
sion, and / « 0.46 in the limit dominated by the bending 
energy (it is the latter limit that is relevant at very small 
H for finite £q). Then, for two hard disks, PFA gives 



PFA 



-1/2 



4/2 



0(x) 



H 



a 

05) 

We plot the ratio of the Casimir energy, Tc [Eq. (|1^)1 
calculated for different cutoff m in the regime dominated 
by the bending energy, and J-pfa [Eq. (|l5|)1 in the inset of 
Fig. 0, and find that the ratio Tq) 'J^fa approaches one 
only at short separations. As the figure shows the higher 
order multipoles are necessary at shorter separations. 

Analytical results for the Casimir free energy ( [T2] ) can 
also be obtained in the weak coupling regime (small A), 
regardless of the separation distance R, or for any A if 
R S> a. These regimes correspond to keeping the first 
term in the expansion of the logarithm in Eq. ([l2]) . For 
the weak coupling regime (small A), we can further sim- 
plify the calculations by replacing G^ with the bare GF, 
Go- The full analytical expressions pll are too cumber- 
some to quote here, and we only present simplified results 
for three important parameter ranges. 

(a) £q ^> R, regime dominated by the bending rigidi- 
ties of the membrane. The Casimir energy ([[2]) has the 
asymptotic form (Eh with n — 4 and the coefficient 



A = (4Bf 



4)Bl 



(16) 



B 9 



K 



4k + kq 



A f 



4(k - kq) + 2(k — k ) 



2k 



K — Kq 



It is remarkable how the Casimir energy depends on the 
flexibility of disks in this regime. For inclusions with 
finite rigidities, Eq. ( |l^ ) is proportional to k — kq = Xk%, 
i.e., a discontinuity in k is required for a non-zero Casimir 
force. The general expression (|l^) reproduces the results 
obtained previously in two limiting cases. In the rigid- 
disk limit [Q, H , where both k and — k are infinite, Eq. (|l|) 
gives A — 6 (horizontal portion of the line shown with 
red symbols in Fig. |l|). In the weak-coupling limit |J, 
it gives A = — A 2 kiKi/2kq. 



4 



(b) 4 < i « fi, regime dominated by the surface 
tension of the membrane. We find that the leading-order 
power law term in a/ R (resulting from dipolc-likc fluctu- 
ations around the inclusions) is zero, and the next-order 
terms give the Casimir energy (|l2|) falling off much faster, 
Tc oc l/R 8 , with the coefficient proportional to The 
full result being too bulky, we only present the strong- 
coupling limit, /3AJ-"c = — 9(a/i?) 8 (in agreement with 
Refs. H and ||), and the leading-order contribution in A 
at weak-coupling, 



PJ$> = 



36A 



36(k — kq) a 



(17) 



Note that we also obtain the same power law in the pres- 
ence of line tension energy on the inclusions boundary, 
in which case Eq. (||) has a ^ ctq inside inclusions Jf9|. 
While a power law oc f / R 8 has been previously obtained 
|^|, ^| for hard inclusions in soap films (kq — ~ko = 0), 
we find it remarkable that the Casimir energy ([H]) de- 
pends on the difference between the Gaussian rigidity of 
inclusions and that of the membrane. 

(c) a < £q < R, with both surface tension and bending 
rigidities of the membrane relevant. The full analyti- 
cal expression |I7J] for the Casimir energy between two 
disk-like inclusions contains terms decaying like an in- 
verse power of the distance R, and exponentially decaying 
terms oc K m (aoR). In particular, for £q <§; R, the terms 
oc K m (aoR) are exponentially small. In this case the 
Casimir interaction energy scales oc l/R 8 , with the coeffi- 
which is a complicated function of parameters, ospo — 



In conclusion, we have developed an exact method 
for computing the Casimir energy between elastic inclu- 
sions of arbitrary shapes embedded in a biological mem- 
brane under tension, characterized by the surface ten- 
sion Co and bending and Gaussian rigidities, kq and Kq. 
The method allows to calculate the Casimir forces in 
all ranges of parameters and for all separations. The 
Casimir energies are fully characterized by the objects' 
"scattering" matrices, which encode the shapes and me- 
chanical properties. In particular, for two elastic disks, 
the Casimir energy scales as oc l/R 4 for R < £q, and 
crosses over to cx l/R 8 for R > £q. At short distances, 
the Casimir energy is large; for hard disks our findings 
agree with the corresponding PFA results, Tq cx iJ~ 1//2 . 
One interesting result is that the Casimir energy is 
strongly suppressed for inclusions whose Gaussian rigid- 
ity k equals that of the membrane. 
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ciall y in the region i 



The exponentially small terms 



became relevant when £ ~ R, where Casimir energy 



J3 



crosses over Lu the snrall-<j regime (a) with Fc oc l/R 



This crossover can be seen in Fig. I . A representative 
case corresponds to In = 10a, where A = — (R/a^J-r, is 



nca jly constant for R < In, is strongly reduced for larger 



R, and eventually crosses over to oc l/R 4 (J^q oc l/R 8 ) 
for R ^> £q . — At smaller R, the same as} r mptotic power 
law lis also seen, e.g., fui £q/u — 1. 



Nble Ihal the distance dependence of Lhe Casimir en- r^r 

erg y is the same, oc l/R 8 , as long as In -C R, which in- [g] 



clu qcs regimes (b) and (c). In the regime (b), dominated 



by the surface tension, this power law can be obtained 
by treating the inclusions as point-like objects in the ef- 
fecti ve field theory (EFT) p. For inclusions that arc free 

ta 



to tf lt with the membrane, the expansion starts with the 
quadrupole terms [||. In the regime (c) the higher-order 
multipolc terms in the EFT expansion diverge as increas- 
ing powers of lp/a ^> 1. However, the contributions to 
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the Casimir energy coming from higher-order multipolc [13 



terras also get suppressed as increasing powers of 1 / R 



As a result, the leading-order quadrupole terms domi- 
nate, which again gives Tc oc l/R 8 for £q <C R. For 
£q ^| h, where we recover the exponentially small terms 
oc it, 

EFT approach cannot be used. 
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i{R/£o), all multipoles contribute equally and the 
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